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We use the iterative unitary matrix multiply method to calculate the long time behavior of the 
resonant quantum kicked rotator with a large denominator. The derealization time is exponentially 
large. The quantum wave delocalizes through degenerate states. At last we construct a nonresonant 
quantum kicked rotator with derealization. 
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Introduction. — The quantum kicked rotator (QKR) [1], 
which describes a periodically kicked rotator, is one of the 
most studied model of quantum chaos [2] . The classical 
correspondence of QKR is the standard map [3, 4]. Clas- 
sically, the energy of the rotator grows without a limit. 
But to a quantum rotator, if the kick frequency and the 
rotator frequency is commensurate, QKR delocalizes in 
the momentum space and if incommensurate, QKR gen- 
erally localizes [1]. Fishman et al explained the classical 
and quantum difference by transforming QKR into an 
Anderson localization problem [5]. 

In the paper, we try to understand how the derealiza- 
tion of commensurate cases happens. This is important 
for several reasons. First, the commensurate (incommen- 
surate) case is described by a rational (irrational) num- 
ber. The qualitative statement that derealization hap- 
pens to the commensurate cases is correct but incom- 
plete. We want to gain a quantitative understanding. 
Second, no physical quantity is rational or irrational. A 
physical quantity has only several significant digits, while 
the distinction between rational and irrational numbers 
depends on infinite significant digits. An infinitesimal 
error can change a rational (irrational) number into an 
irrational (rational) number. While we expect the system 
changes little from our experiences of studying physics as 
was emphasized by Hofstadter [6] . To recognize and rec- 
oncile the conflict is one aim of quantum chaos. Third, 
Fishman et aVs result [5] seems to tell us localization 
happens to all the incommensurate cases. Is there at 
least one incommensurate case for which derealization 
happens? Casati et al has derived a quantum Lyapunov 
equation to describe the difference between the dynamics 
of commensurate and incommensurate cases [7, 8] . Based 
on the formula, Casati et al claimed there are some in- 
commensurate cases of derealization [7]. But their ar- 
gument is problematic [8] from the perspective of the 
exponentially large derealization time discovered in the 
paper. 

In the paper, we prove by numerical calculation for 
the commensurate case with a large denominator, the 
derealization time is exponentially large. Such a large 
denominator effect is explained by the degenerate per- 
turbation theory, which is based on the observation that 
degenerate states are the derealization path. Localiza- 
tion of incommensurate cases can be understood to be 



caused by the large denominator effect. The large de- 
nominator effect and the quantum Lyapunov equation 
[7, 8] partially reconcile the conflict between commen- 
surate and incommensurate cases and naturally lead to 
an incommensurate case of derealization. This partially 
solves the problem: to find an incommensurate case of 
derealization, posed by Casati et al [7, 9] and gives a 
counterexample to Fishman's argument [5], although a 
very weak one. 

Numerical methods. — For a system with a periodical 
Hamiltonian, the unitary operator of one period is the 
Floquet operator F. The unitary operator of 2 N periods 
is F 2 . 
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From Eq. (1), we can calculate F 2 from F by iter- 
atively multiplying the unitary matrices for N times. 
This method is referred as the iterative unitary matrix 
multiply method (IUMM). It is impossible to calculate 
very long time behavior of QKR using the usual fast 
Fourier transform method [5]. IUMM is actually the 
same method as direct diagonalization or the matrix vec- 
tor multiply method used in the original paper [1] of 
QKR. See the section III and IV of [1]. 

Calculation results. — The Hamiltonian of QKR is 
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where h is the Planck constant, r the kick period and k 
the kick strength. The matrix element of F is 



(n|-F|m) = exp(— ifir- 
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where \n) = ^=e me . We apply IUMM to QKR. In the 
calculation h = 1, k = 1 and r = ^ = ^g. The initial 
state is |0). This is the commensurate/resonant case with 
a large denominator q = 20. F nm 
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FIG. 1: QKR wave function at different time. TV is at time 2 n t. n is \n) and c n is the base-10 logarithm of 
the absolute value of the wave function on \n). n is from —500 to 500 in our calculation. 




FIG. 2: Clearer figures of QKR wave function. Note the changes of peaks and valleys from N — 26 to TV — 43. 
The peaks indicate what \n)s actively contribute to the wave propagation. 
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every n and m. The rotator will delocalize in the future, 
nevertheless it delocalizes very slowly. 

FIG. 1 and FIG. 2 show the same calculation results. 
The distribution of QKR wave function is clearer in FIG. 
2. Before N = 20 (the time is 2 20 t = 1.05 x 10 6 r), 
QKR docs not delocalize at all. To a resonant case, this 
is unexpected. At N = 20, peaks (local maxima) and 
valleys (local minima) appear. Naively one expects, from 
iV = 20 to 50, peaks should be at \n = 20 x 7)s, where 7 
is an integer, because such states are resonant with the 
|0). Actually \n = 20 x 7)s are always valleys. Before 
N = 39, the peaks are at \n = 10 x O ± 3)s, where O 
is an odd integer, such as n — 13, 27, 33, 47, 53, • • • . At 
N = 50, the peaks are at \n — 20 x I ± 3)s, such as 
n= 17,23,37,43, 

At N = 26, the triangle like wave function outside 
n = forms and at N = 30 it flattens. The height of 
the flattened wave function is approximately 10 -4 . The 
main distribution is still at |0). From N = 30 to 39, 
the needle like wave function around n = becomes a 
triangle like one. From N = 39 the triangle expands and 
at N = 45 the wave function flattens again. The height 
of the QKR wave function is approximately 10~ 2 . At 
N = 50, a triangle like wave function forms again. 

Around n = 0, the circulation: needle — ► triangle — > 
expanded triangle — ► flattened, drives the whole delocal- 
ization process. Will the wave function around n = be 
totally flattened by one more circulation or several more 
after N = 50? Will \n — 20 x 7)s finally become the only 
peaks when N is very large and how? We do not know! 
Our calculation overflows around N = 60, which may 
be caused by the ununitarity of the truncated Floquet 
operator. 

At what N should the wave function be considered as 
delocalized? At N = 30, there is still lots of distribu- 
tion of the wave function at n = 0. N = 45 is more 
proper than N = 30. We define T T / 27r x r as the delo- 
calization time of QKR with the period r. For simplic- 
ity, we also refer T T / 27T as the delocalization time. For 
t = 2ir/q, QKR delocalizes exponentially slowly. We es- 
timate T\i q ~ exp(cg/fc), where c is a factor that depends 
weakly on k and q. If we assume QKR is delocalized at 
N = 45, e 20c r w 2 45 t and c w 1.56. 

Delocalization path and degenerate perturbation the- 
ory. — The delocalization time can be estimated from 
the degenerate perturbation theory. The sequence 
{— ^n 2 Mod(l)}„ = o,i,... .20, which are the phases of F nn 
divided by 2n, is 
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In a period, there is four ||s, four |^s, two y^s, two |s, 
two |s, two |s, two t^s, one \ and one 0. The quantum 
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wave is easier to propagate between |is or between 
So |n = 20 x I)s are valleys from N — 20 to 50. When 
n = 1,9,11,19,21,29,31,39, the phases are |§ x 2tt. 



The intervals between two degenerate states are 1 and 
8. When n = 3,7,13,17,23,27,33,37, the phases are 
|i x 2n. The intervals are 4 and 6. So the wave is easiest 
to propagate between n = 3, 7, 13, 17, 23, 27, 33, 37, which 
are peaks in FIG. 2. But we do not know why peaks are 
at only some of the |n = 10x/±3)s. Peaks even change 
from the \n = 10 x 0±3)s to the |n = 20 x J±3)s as we 
discussed above. 

From the degenerate perturbation theory, if the wave 
propagates through the path |0) — ► |20), F is approxi- 
mated by 
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The eigenvalues of F appr is J (l) ± J2o(l)- So after ap- 
proximate Jo(l)/^20 = 1-98 x 10 24 -time kicks, the wave 
function will be transferred from |0) to 1 20) . This is far 
larger than 2 45 = 3.52 x 10 13 of our numerical result. A 
more exact estimate has to take into account other degen- 
erate states. The quantum wave can propagate through 
the path |0) |1) -» |9) -> |11) -» 1 19) -» |20). The 
states contributing to the wave propagation are mainly 
these states. So F is approximated by T'appr, which only 
considers the states in the delocalization path. 
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The propagation time from |0) to |1) is J (l)/Ji(l); from 
|1) to |9> is J (l)/J 8 (l); from |9) to |11) is J (l)/J 2 (l); 
and so on. The delocalization time from |0) to 1 20) is 
estimated to be 



Ti/ q 



Jo(l)Jo(l)Jo(l)Jo(l)Jo(l) 
Ji(l)J 8 (l)J2(l)J 8 (l)7i(l) 



= 1.33 x 10 15 . (7) 



This is more realistic than Eq. 5. Another path of wave 
propagation, |0) -» |3) -» |7) -» 1 13) -» |17) -» |20), 
gives 



Mi) Mi) Mi) Mi) Mi) 

J 3 (l)J 4 (l)J6(l)J4(l)^3(l) 



5.34 x 10 1 
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which is close to the delocalization time 2 45 = 3.52 x 10 13 . 
One problem of the degenerate perturbation theory is 
F appr is not unitary. 

An incommensurate case of delocalization. — The 
smaller q, the faster the delocalization. If t/2tt = p/q w 
p'/q' and q' <C q, QKR with r = 2irp/q delocalizes 
quicker because it is closer to a stronger resonance. But 
r = 2ir/q is far from any strong resonance in all the 
(r = 2np/q)s, where p = 1, 2, • • • , q — 1, q. So it has the 
largest delocalization time and 7\/ 9 w exp(cg/fc) is the 
upper limit of delocalization time in all the (r = 2irp/q)s. 
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Now we construct irrational t/2tt with de-localization. 
Imagine two QKRs with almost equal kick period t and 
t' and the equal kick strength k. St = |r — r'| <C 1. 
U(M,t) is the M-period unitary operator with the kick 
periods r and U(M,t') with r'. The difference between 
the matrix elements of two unitary operators [7, 8] 

\U(M, T ) nm - U(M, T ') nm \ < 1 M*k 2 6T. (9) 

As the particular value of 7 is not important, we set 
7 = 1. Before {e/{k 2 St)) 1 /* -time kicks, \U{M,r) nm - 
U(M,r') nm \ < e. 

We consider k = 1 and construct 

r/27r =l/q+ l/[exp(3cig)J + l/[cxp(3c 2 exp(3cig))J 
+ l/[exp(3c 3 exp(3c 2 exp(3cig)))J 
+ ••• 

+ l/[exp(3c„ • • • exp(3c 3 exp(3c 2 exp(3cig))))J 
+ ••• . 

(10) 

[xj is an integer around the real number x (For the con- 
venience of the argument below, \x\ is not the same as 
the floor function in mathematics.) and ensures every 
term is a rational number, q is a positive integer such as 
20. ci > cid and cid is the factor in the derealization 
time T\/ q — exp(cidg). c 2 > c 2 d and c 2 d is the fac- 
tor in the derealization time T 1 / qi = exp(c 2d gi), where 
qi = [exp(3cig)J. And so on. 

|T/27r-l/g|wl/exp(3cig). (11) 

After ^/Tff-timc kicks, the dynamics of r and r' = 2n/q 
will not diverge from each other much due to Eqs. (9) 
and (11). QKR with r propagates to a domain in the 
momentum space as large as l\. We choose c\ ^> Cid to 
ensure l\ ^> fc 2 /4 = 1/4. We choose [exp(3ci(7)J to be an 
integer approximately exp(3cig) and to be multiples of q. 
So from ^fqi kicks to ^fqi kicks, the derealization speed 
of QKR is larger than or equal to QKR with 2-k jq\. After 
^/g^-time kicks, QKR propagates to a larger domain I2 > 
l\. We choose c 2 S> c 2 d to ensure Z 2 3> h- And so on. 
= 00. So QKR with the kick period r will delocalize. 



Even if P 1 / q ^ exp(cq/k), we can always construct 

QKR with t in Eq. (12) delocalizes. 

We note similar irrational numbers have been con- 
structed by Avron et al concerning the Harper equation 
[10] and by Berry [11] and Prange et al [12] concerning 
the Maryland model. It cannot be a coincidence that 
similar numbers are constructed to three totally different 
problems. We think such irrational numbers universally 
have similar behavior with rational numbers in problems 
of quantum chaos. The way to construct irrational num- 
bers in Eqs. (10) and (12) is very general and our argu- 
ment depends on Eq. 9, which is a universal quantum 
Lyaponov equation [8] . 

Problems. — Some problems remain. First, how does k 
influence Tjy 9 ? Second, how to estimate T p / q ? Do T p/q s 
generally approximate to T\/ q l Third, is there one in- 
commensurate case of derealization, which is not similar 
to Eqs. (10) and (12)? We think localization happens to 
the general Liouville number t/2tt, such as the Liouville 
constant. 

Conclusion. — First, we have calculated the long time 
behavior of QKR using IUMM. It is discovered the der- 
ealization time is exponentially large for large denomina- 
tors. Second, we have constructed an irrational number 
of derealization. Concerning QKR, Eqs. (10) and (12) 
are the first irrational number with derealization ever 
known. Both results have important meaning for the 
theory of QKR. Third, the large derealization time is 
explained by the degenerate perturbation theory, which 
is suggested by and consistent with the derealization 
path of the numerical calculation. The phenomena that 
the wave propagates between degenerate or almost de- 
generate states may be found in many other systems. 

This work is supported by the National Natural 
Science Foundation of China under Grant Numbers 
10674125 and 10475070. I would like to thank Profes- 
sor Fishman for helpful discussions. 



[1] G. Casati, B. V. Chirikov, F. M. Izraelev, and J. Ford, in 
Stochastic Behavior in Classical and Quantum Hamilto- 
nian Systems, edited by G. Casati and J. Ford (Springer, 
Berlin, 1979), vol. 93 of Lecture Notes in Physics, pp. 
334-352. 

[2] H. J. Stockmann, Quantum Chaos: An Introduc- 
tion (Cambridge University Press, Cambridge, England, 
1999). 

[3] B. V. Chirikov, Phys. Rep. 52, 263 (1979). 

[4] J. M. Greene, J. Math. Phys. 20, 1183 (1979). 

[5] S. Fishman, D. R. Grempel, and R. E. Prange, Phys. 

Rev. Lett. 49, 509 (1982). 
[6] D. R. Hofstadter, Phys. Rev. B 14, 2239 (1976). 



[7] G. Casati and I. Guarneri, Commun. Math. Phys. 95, 
121 (1984). 

[8] T. Ma, General theory of the quantum kicked rotator. I, 

nlin/0710.1661. 
[9] G. Casati, J. Ford, I. Guarneri, and F. Vivaldi, Phys. 

Rev. A 34, 1413 (1986). 
[10] J. Avron and B. Simon, Bull. Am. Math. Soc. 6, 81 

(1982). 

[11] M. V. Berry, Physica 10D, 369 (1984). 
[12] R. E. Prange, D. R. Grempel, and S. Fishman, Phys. 
Rev. B 29, 6500 (1984). 



